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INFERENCE OF THE BOUNDARY LAYER STRUCTURE OVER " ‘ 

THE OCEANS FROM SATELLITE INFRARED MEASUREMENTS 

INTRODUCTION 

Temperature profile in the atmosphere has been successfully derived from passive infrared 
and microwave radiation measurements made from remote platforms (see for ex. Wark and 
Fleming, 1966; Conrath et aL, 1970; Staelin et aL, 1975; Smith and Woolf, 1976). In the 
infrared, spectral measurements in the vibration rotation bands of CO 2 at 1 5 and 4.3 /un are 
used for this purpose, while in the microwave the measurements in oxygen, O 2 , bands around 
55 GHz (0.5 cm) are used. Both these gases. CO 2 and O 2 , are uniformly mixed in the atmo- 
sphere and htnce their concentration decreases, with height, with a scale height of about 8 km. 
As such the weighting function for different spectral regions in the absorption bands of CO 2 
and O 2 tend to be broad. For this reason the vertical resolution of the remotely sensed tem- 
perature profile from CO 2 and O 2 bands is rather poor. In particular the vertical resolution 
of the temperature profile in the lower troposphere is about 5 km which results in smoothing 
out of some meteorologically important small scale temperature features. : 

Water vapor in the atmosphere has a uch smaller scale height of about 2 km and is heav- 
ily concentrated in the lowest layers of the atmosphere. However, unlike CO 2 and O 2 , ver- 
tical distribution of water vapor in the atmosphere is highly variable, particularly over land, 
which makes it difficult to use water vapor as an optically active gas to do remote sensing of 
temperature. This difficulty is substantially removed on large water bodies where it is possi- 
ble to model the mean relative humidity profile with the help of climatology. 

This mean relative humidity profile is significantly modified in the intcrtropical conver- 
gence zone (1TCZ) and in the trade wind inversion regime. When a low level convergence is 
present in the atmosphere, a deep convective layer is formed and this layer tends to have a 
high relative humidity. On the other hand if stable conditions in the temperature profile, _ 
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such as an inversion, exist, the water vapor is not convected up freely, and above the inver- 
sion usually a dry layer with low relative humidity exists. In this fashion strong correlation 
between the temperature and water vapor profiles is developed over water bodies depending 
on the conditions prevailing in the boundary layer. Given this a priori information of the 
correlation we can infer from the water vapor spectral measurements some details in the tem- 
perature profile close to the surface. In this study we have followed this approach to remotely 
sense the characteristics of the boundary layer of the atmosphere over the global oceans be- 
tween about 50°N to 40°S for three different periods, about 3 months each, during the year 
1970. The spectral measurements made by Nimbus 4 Infrared Interferometer Spectrometer 
(IRIS) are used for this purpose. 

EXAMINATION OF THE INFORMATION CONTENT IN IRIS SPECTRA 
The Nimbus 4 IRIS gathered spectral measurements from 400 cm* 1 to about 1400 cm* 1 
over the globe for a period of about one year (April 1970 to Jan. 1971). The spectral reso- ~ 
lution of this instrument is 2.8 cm* 1 and the noise in the spectral data is about 0.5 erg cm* 1 
sr* 1 s* 1 (Hanel et al., 1972). The field of view of the IRIS is about 95 km in diameter. In 
Fig. 1 a typical brightness temperature spectrum of IRIS is shown. The 15 pm CO2 band 
and 9.6 pm O3 band are distinctly shown in the spectrum. The various water vapor bands, 
the rotation band around 20 pm, the window regions at 1 1 pm and 9 pm, and the 6.3 pm 
vibration rotation band, show some spectral details at 2.8 cm* 1 resolution. These spectral 
details are produced by the water vapor lines which are distributed in some random fashion. 
Whenever a few strong water vapor lines come close together a strong minimum in the bright- 
ness temperature can be noticed. This implies that even at 2.8 cir* 1 resolution we have some 
crude information of the water vapor lines retained in the spectiai data. On the other hand 
the line information in the CO2 and ozone bands is almost completely smeared out at 2.8 cm" 1 

resolution as the spectral lines in these bands are more closely arranged compared to the 

water vapor bands. 
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Figure 1. Nimbus 4 IRIS brightnes emperature spectrum taken over the equatorial Pacific 
(2.7°N, 140.1°W, April 20, 1970). Notice the crude structure of the water vapor lines at 
2.8 cm* 1 resolution in the 9/im region (1100 - 1230cm* 1 ) 






The sharp maxima (peaks) aiid minima (valleys) in the water vapor bands suggest that the 
absorption coefficient of water vapor changes significantly in narrow spectral intervals. This 
implies further that radiances, corresponding to the brightness temperature, at the peaks and 
valleys in the spectrum, arise from significantly different altitudes in the atmosphere. 

In this study we are interested in examining the information pertaining to the lowest layers 
in the atmosphere. Hence we have considered the water vapor spectral information in the win- 
dow regions of 1 1 -1 3 pm and 8-9 pm. These window regions are referred to as 1 1 pm and 
9 pm regions respectively. In the 9 pm region the continuum absorption due to self-broad- 
ening effect is weaker (Bignell. 1970; Burch, 1970) while the lines are stronger as compared 
to the 1 1 pm region. For this reason we have considered the 9 pm region for the purpose of 
the present study. In order to appreciate the nature of the water vapor lines in the 9 pm 
region we have shown in Fig. 2a the transmission function of the water vapor, for 2 gms of 

preripitable water in the atmosphere, produced by a detailed line by line calculation at res- 

olution of 0.1 cm" 1 (Kunde and Maguire, 1974). Fig. 2b shows the transmission function 
when degraded to a resolution of 2.8 cm" 1 . revealing that the water vapor line information 
is not obliterated at this resolution. 

From an examination of Fig. 2b and the 9 pm window region shown in Fig. 1, we infer 
that in the peaks, produced by weak continuum absorptions; information close to the sur- 
face is being sensed. This corresponds, to a close approximation, to the surface temperature. 
The valleys on the other hand, formed by moderately strong lines, sense an effective temper- 
ature some what higher up in the atmosphere. The difference between the peaks and valleys, 
which ranges from about 1 to i2°C in temperature, as will be shown later, gives us some in- 
formation about the structure of the atmosphere in the first 2 or 3 km near the surface. 



WAVE LENGTH ( M m) 


Figure 2b. Same as Fig. 2a when degraded to 2.8 cm*' resolution. 
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THEORETICAL CONSIDERATIONS ““ ‘ 

It was pointed out in the previous section that the water vapor spectral measurements in 
the 9/im window region could yield information about the atmospheric stratification in the 
boundary layer over theoceans. We will develop this idea using the radiative transfer .ry. 

The intensity I„ measured at the top of the atmosphere, at wave number v y in a non- 
scattering atmosphere in local thermodynamic equilibrium, is given by 

I* = B„(p 0 ) r„(p 0 ) + f B„(p) d r„(p) (1) 

where 

B is the Planck intensity, 

T(p) is the transmission of the atmosphere from any given pressure p to the top of the 
atmosphere, 

p 0 is the surface pressure. 

When we apply this equation to the spectral data in the 9 /im water vapor window region, 
the first term on the right hand side, i.e. the contribution from surface is essentially governed 
by the surface temperature and the total precipitable water content in the atmosphere. The 
second term accounts for the emission of the whole atmosphere, and it depends on the man- 
ner in which the water vapor and temperature are distributed. Within the 9 /am window, as 
shown in the IRIS spectrum (Fig. 1), at the peaks where the absorption is weak the contri- 
bution from the second term is much smaller than the first one. Thus the intensity at these 
peaks depends only weakly on the atmospheric stratification. However, in the valleys of the 
spectrum, where the absorption is considerably larger, the second term in the equation ex- 
ceeds the first one. Now if we take the difference A1 between the intensity at the peak. Ip, 
and the intensity I v in the valley, we get some information about the atmospheric structure. 




An equation for AI may be derived fronT equation (1) as follows' 

/ o^ 

^(p)d[Tp(p)-r v (p)] (2) 

l T p(Po)- r v(Po)] 

where 15 is the average Planck intensity corresponding to the peak and valley in the spectrum 
which are spaced within 10 cm" 1 

and T p and r v are the transmissions in the peak and valley of the spectrum. 

From equation (2) we see that AI will be zero when the atmosphere is isothermal having a 
temperature equal to the surface temperature. On the other hand if the temperature increases 
from the surface to considerable height AI can be negative. Such negative values are com- 
monly noticed in the polar latitudes. 

Equation (2) may be generalized as 

£ = f[T s ,T(p),R.H.(p)] (3). 

where the line strength £ is the brightness temperature difference corresponding to AI, , 

i 

T(p) and R.H.(p) are the temperature and relative humidity distributions as a function of j 
pressure, p, in the atmosphere, 

T s is the surface temperature. 

Now if we are given the mean conditions of temperature profile, T(p), associated with a given 
surface temperature T s and relative humidity profile, R.H.(p), we can express the line strength 
£ corresponding to these mean condition as 

5 = f [T s ,T(p),CT.(p>] (4) 

When the temperature and the relative humidity profiles in the atmosphere for a given T s 
differ from such mean conditions the line strength £ departs from ? by an amount, say A£. 
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Such a perturbation A8 can than be expressed as 

A8 = 8 - 8 = g [ T s , T'(p), R.H.'(p)] (5) 

where T' and R.H.' are deviations from average. 

The surface temperature T s can be determined independently from the 1 1 pm window 
measurements of IRIS using a method develop*, by Prabhakara et al., (1974). Thus from 
equation (5) we see that from o.ie quantity, A8, we can infer one dominant mode or property 
of the atmosphere implicitly containeu in the right hand side of the equation. On the water 
bodies, as will be discussed later, the intimate relationship between the characteristics of the 
temperature and water ,f apor profile in the boundary layer can be modelled to define such a 
dominant mode. 

SPECTRAL SIMULATIONS FOR DIFFrRF>IT ATMOSPHE^’C MODELS 

_ In order to study the manner in which 8, as defined previoi sly, changes as a function of 

the atmospheric conditions we have developed a radiation computational scheme that can 
simulate the spectral d^.a with a resolution of 2.8 cm' 1 in the water vapor absorption regions 
from 400 cm -1 to 1400 cir ' The transmission function, r, of the water vapor needed in 
these Iculations is taken to be the product of three components: rj associated with water 
vapor lines, Tp produced by the continuum due to foreign broadening, and r e . suiting from 
the e-type absorption (Bigncll, 1970). The transmission function of the water vapor lines 
T{ is derived using a multiple regression scheme similar to the one proposed by Smith (1969). 
The details of the foreign broadening and the e-type given by Kurde and Maguire (1974) arc 
adopted. 

Th:s simulation program is used to synthesize spectral data for several different atmospheres 
ranging from tropics to high latitudes. The atmospheric temperature and relative humidity 
distriDutions corresponding to mean conditions in the latitudinal belts 0 to 10°, 10 to 30°. 
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and 30° to 50° on both the hemispheres for four seasons are derived by Nieman (1977), ~ 
from radiosonde data measured during the year 1972, over the globe (see Table 1 ). Some of 
these mean profiles of the temperature and relative humidity are shown in Fig. 3. It is inter- 
esting to note that in all cases the mean relative humidity profile is very similar and decreases 
from surface up to about 200 mb almost linearly as a function of height. These data are 
heavily weighted by land stations. However, from a limited sample of ship and island radio- 
sonde data taken over north Atlantic (see Tables 6a and 6b) during the year 1970, we find a 
mean value of 81% relative humidity near the surface, and a similar decrease with height. For 
this reason we have adopted for the radiative transfer simulations, one relative humidity pro- 
file. as shown in Fig. 3, for all oceanic regions. 

From the synthetic spectral data, calculated with the simulation program, using the above 
mean atmospheric conditions, we have obtained the brightness temperature at eight peaks in 

the 9/im region - 1 127.7, i 140.2, 1 158.3, 1 168.0, 1 181.9, 1 193.1, 1202.8, 1232.0cm" 1 - 

and averaged them. Similarly the brightness temperature at eight adjacent valleys - 1 136.1, 
1149.9, 1165.3. 1173.6, 1186.1, 1197.2, 121 1.2, 1225.1cm" 1 - are averaged. The differ- 
ence between these two averages gives an estimate of the line strength representative of the 
spectral region 1 1 25 to 1 235 cm" 1 . This value of £, in brightness temperature, is shown in 
Tuble 2 for the various mean atmospheres. In Fig. 4 we show graphically the relationship 

i 

between the surface temperature of the various model atmospheres and ?. The total precip- ■ 
table water content w is also shown in this figure as a function of the surface temperature. 
Both w and C, as seen from the figure 4, increase with T s in an analogous fashion. Tire in- ; 
crease in w may be readily understood in terms of the moisture holding capacity of the 
warmer •Umosphercs. Tire growth of ? needs some explanation. i 

The behavior of C can be explained with the help of the weighting functions, dr/dz, asso- 
ciated with the valleys and peaks in the spectra. In Fig. 5 the weighting function corresponding 
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/ Figure 3. Mean relative humidity and temperature profiles for tropics and midlatitudes summer and winter 
(from Nieman, 1977). The adopted R.H. profile is also shown. 
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Table 2 


Mean line strength ? and total precipitable water w as a function of the 
surface temperature T s for different mean atmospheres. 



t s(K) 

I(°C) 

w (g/cm 3 ) 

Tropics 
10°N - 10°S 

298.9 

9.4 

4.28 

Subtropics Summer 
10°N - 30° N 

298.6 

9.3 

4.29 

Subtropics Winter 
10°N-30°N 

291.9 

6.6 

2.67 

i 

Temperate Summer 
30° N - 50° N 

292.8 

7.1 

2.91 

Temperate Winter 
30° N - 50° N 

273.9 

2.8 

0.92 


















LINE STRENGTH <°C) 



SURFACE TEMPERATURE (K) 

Figure 4. Dependence of the 9 pm mean line strength V and the total water vapor content w on the surface 

temperature T s over the oceans. 
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WATER VAPOR CONTENT (g/cm 2 ) 






to a valley at 1225 cm -1 and adjacent peak at 1232 cm -1 are shown for three atmospheres. 
The weighting function associated with the valley has a distinct maximum which moves 
down from about 3 km to 1 .5 km as the total water content decreases. The weighting func- 
tion at the peak on the other hand does not show such a feature. 

The height at which the maximum in the weighting function occurs is around unit optical 
depth in the atmosphere. At 1225 cm' 1 , a valley in the spectrum, unit optical depth is reached 
at a higher elevation in a warmer and, therefore from climatology, more humid atmosphere. 

In the peak at 1232 cm -1 unit optical depth is not reached in all the cases. 

The strength of the line 5 is thus crucially dependent on the super incumbent water vapor 
in the troposphere. This relationship between 8 and w is shown in Fig. 6 by a solid line. 

This pseudo linear relationship suggests that when mean conditions of temperature and hu- 
midity profiles are present £ grows in proportion to w. 

From the above discussion we may conclude that, over the oceans, the relationships shown. 

in Figs. 4 and 6 represent basic correlations among the surface temperature T$, w and 8. 


These correlations may be broken when atmospheric temperature and water vapor profiles : 

i 

are perturbed from mean conditions. Such perturbations are present over the global oceans. 

i 

In the following discussion these climatologically significant perturbations are examined ' 

i 

in some detail. i 

Trade wind inversion J 

— 1 

Trade wind inversion is a climatological feature associated with oceanic subtropical anti- 


cyclones. In the northern hemisphere the axis of these anticyclones is tilting from NE toSW 


(see for ex. Haurwitz and Austin, 1942) and thereby the vertical motion associated with 


these anticyclones has an organized pattern. On the NE side there is a pronounced subsidence 
while on the SW the vertical motion is weak. This organized motion of the subtropical anti- 


cyclones produces strong inversion on the NE side. The observations of the Meteor oceanographic 

i 




0 


2 


4 


6 


WATER VAPOR CONTENT (g/cm 2 ) 

Figure 6. The relationship between the 9pm line strength and the total precipitable water for mean atmo- 
spheric conditions, for inversion conditions (from Tables 4a, b and c), and ITCZ cases (see Tables 5a and b). 
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expedition (Ficker. 1936; Riehl, 1954) reveal the general nature of the trade wind inversion 
over the North and South Atlantic. The more recent experiments related to GATE, in par- 
ticular BOMEX and ATEX (Augstein et al., 1973; Dunckel et al., 1974), reinforce the obser- 
vations of the Meteor expedition. 

In Fig. 7 an example of trade wind inversion is presented to illustrate the height of the in- 
version. the temperature increase AT from bottom to top of inversion and the corresponding 
decrease in relative humidity AR.H. Ficker (1936) has constructed, from Meteor data, maps 
of the change in temperature and relative humidity from bottom to top of the inversion, as 
well as a map of the height of the inversion. Riehl (1954) shows these maps and presents a 
lucid discussion. A close examination of thest maps reveals intimate connection between the 
height of the inversion, the temperature change, and relative humidity change associated with 
the inversion. In Fig. 8 we show the manner in which the inr-'ase in temperature, AT, and 
the decrease in the relative humidity, AR.H., across the inversion, are related. This shows a 
good positive correlation between the two variables. In Fig. 9 the AT and the height of the 
inversion are related. This relationship is not as well defined as the previous one although a 
negative correlation between the two variables is suggested. These correlations are developed 
in a consistent fashion by subsidence motion. That is to say the stronger the subsidence the 
lower the height of the inversion and the larger the temperature increase from bottom to top 
due to adiabatic warming (Riehl, 1954). The change in relative humidity is primarily due to 
the capping of convection by the inversion, which limits the height of the lower humid layer, 
while above the inversion sinking of dry air and dynamic warming lower the relative humid- 
ity. In this manner a correlation between the temperature profile and the water vapor pro- 
file in the troposphere is developed on oceanic regions where the inversion conditions prevail. 

The Indian Ocean expedition has revealed the inversion characteristics over the Arabian 
sea. From the measurements made along 50°E by the ship R.V. Discovery, Ramage (1971) 
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Figure 7. Example of temperature and relative humidity profiles in the trade wind inversion regime 
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AT "C-TEMPERATURE INCREASE ACROSS INVERSION 


Figure 8. Relationship between the temperature increase, AT, and the decrease in relative humidity AR.H 
from bottom to top of inversion, derived from the observation of the Meteor expedition (Ficker, 1936). 
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AT °C-TEMPERATURL INCREASE ACROSS INVERSION 

. figure 9. Relationship between the height h of the inversion and the increase in temperature from bottom 
to top of inversion derived from the observations of the Meteor expedition (Ficker, 1936), 
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shows a meridional cross section of the temperature and relative humidity distribution from' 
7°N to 1 2°N during August 1964. This cross section shows a significant temperature inver- 
sion, up to as much as 10°C warming at 10°N, between the surface and 1 km level. In another 
figure Ramage { 1966) shows, from ship observations made during summers of 1963 and 1964, 
the distribution of the level of inversion over the Arabian sea and the adjoint* . Indiar Ocean. 
The salient features of this distribution are that the height of the inversion is gene -ally low, 
at about 1 km, close to the coast line of Somali, Arabia, Iran and Pakistan. Away from this 
coast line the height of the inversion rises to about 3 km as one goes towar ndian Ocean. 
In a broad sense the observations of the Indian Ocean Expedition reinforce ^cas derr ed 
from the Meteor Expedition. 

The graphical relationships shown in Figs. 8 and 9 enable us to model the inversion charac- 
teristics over the oceans with one parameter say AT which is coupled to AR.H and h. In 

Table 3, three values of the inversion parameter AT and the associated AR.H and h, derived 

from Figs. 8 and 9 are shown. These inversion characteristics are applied to tl e various mean 
atmospheres shown in Table 1. The temperature profile above the inversion layer is obtained 
by a straight u ;e fit between the temperature at 500 mb and the temperature at the top of 

Table 3 


Three parameter representation (Fig. 7) of the trade wind inversion based on 
Meteor data (see Figs. 8 and 9). 


Height of Inversion 

Temp. Increase 

Decrease in Rci. Hum. 

h (km' 

AT (°C) 

AR.H (%) 

1.70 (825 mb) 

4m 

32 

1.15 (880 mb) 

4 

48 

0.75 (930 mb) 

6 

62 


i 
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the inversion. The adopted relative humidity profile shown in Fig. 3 is suitably modified to ' 
represent the inversion conditions. Above the level of inversion the relative humidity profile 
is extrapolated assuming a constant value til) it meets the undisturbed profile. In Fig. 10 these 
three relative humidity profiles of the trade wind inversion a*-- shown, \pplying these three 
inversion characteristics to the tropical, mid latitude sv nmer and the mid latitude winter 
atmospheres (see Table 1 ) we have calculated 5, the strength of the lines in tne 9 jum region. 

In Tables 4a, b and c, the strength of 2 and the total precipitable water w for the three in- 
version c editions are listed. These tables show the dependence of 1! on the boundary layer 
parameters h, AT and AR.H in a joint fashion. We have attempted to decouple the depen- 
dence on these parameters and examine the sensitivity of 2 separately to ti.e temperature 
profile, the water vapor profile and the height of the inversion level. Such an examination 
shows that about 85% of the change in 2 is due to the variation in the water vapor profile, 

while the temperature profile accounts for about 10% and h accounts for about 5%. The 

above analysis shows the overwhelming importance of the water vapor profile. Any inference 
of the temperature profile from 2 thus has to be deduced by way of statistical correlation 
with respect to water vapor profile. 

For a given surface temperature the total water content in each one of these models varies 
as the boundary layer parameters change. The water content is at a maximum when the in- 
version is absent and least when the inversion is at the lowest levels as shown in the Tables. 

As the total water is an integral property reflecting the inversion characteristic we may relate 
it to 2. 

In Fig. 6 the values of 2 and w taken from Tables 4a, b, c, are plotted. It can be seen from 
this figure that all the points are very close to the curve corresponding to mean atmospheric 
conditions. As a consequence of tills result it is possible to infer w from 2 in the presence of 
the inversion condition. 
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PRESSURE (mb) 



RELATIVE HUMIDITY (%) 

. Figure 10. Adopted mean relative humidity profile (shown also in Fig. ?) and the models of the relative 
humidity profile for three conditions (see Table 3) of the inversion. 
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Table 4 


Computed 9 pm line strength corresponding to three trade wind inversions (see Table 3) 
for tropics (a), mid-latitude summer (b) and mid-latitude winter (c). 

(a) Tropical mean atmosphere (T s * 298.9 K, £ = 9.62°C, w = 4.28 g/cm 2 ) 


h (km) 

AT (°C) 

AR.H (%) 

«(°C) 

w (g/cm 2 ) 

(w - w)/w 

1.70 

2 

32 

8.75 

3.71 

0.13 

1.15 

4 

48 

7.30 

2.91 

■H 

mm 

6 

62 

5.22 

1.98 

0.54 


(b) Midlatitude Summer (T s = 292.6 K, 7 = 7.4 8°C, w = 3.35 g/cm 2 ) 


h (km) 

AT (°C) 

AR.H (%) 

£ (°C) 

w (g/cm 2 ) 

(w - w)/W 

1.70 

2 

32 

6.51 

2.89 

0.14 

1.15 

4 

48 

5.18 

2.22 

0.34 

0.75 

6 

62 

3.54 

1.46 

0.56 


(c) Midlatitude Winter (T s = 273.8 K, £ = 3. 10°C, w = 0.96 g/cm 2 ) 


h (km) 

AT (°C) 

AR.H (%) 

£(°C) 

w (g/cm 2 ) 

(w - w)/w 

1.70 

2 

32 

2.62 

0.81 

0.16 

1.15 

4 

48 

1.97 

0.65 

0.35 

0.75 

0 

62 

1.22 

0.43 

0.57 


The surface temperature T s on the ocean can be measured independently from the 1 1 pm 
IRIS window data. Given the surface temperature we can estimate the total water content 
w corresponding to mean atmospheric condition, from Fig. 4. Now if we can combine the 
parameters of the inversion in a model, as discussed earlier, we can relate (w - w)/w to the 
strength of the inversion parameters. We find the relationship applicable to AT is 
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(w - w)/w = A AT 


( 6 > 


where A is a constant approximately equal to 0.1 °C“* . 

This relationship is shown in Fig. 1 1. 

Inter Tropical Convergence Zone. 1TCZ 

Another well known phenomenon over the tropical oceans is the Inter Tropical Conver- 
gence Zone, ITCZ. This is the region where the low level winds converge to produce a rising 
motion which carries water vapor upward and produces above average humidity conditions 
aloft. The temperature inversion above the boundary layer in the ITCZ is generally absent 
or weak (Estoque, 1975). Thus the temperature profile in the ITCZ areas does not show a 
significant departure from mean atmospheric conditions. However the water vapor profile 
reflects the effects of rising motion, resulting in a significant increase in total water content. 

For this reason we can not develop a simple model of the coupling of the temperature and 

water vapor profiles in the ITCZ. 

Observations of the temperature and humidity corresponding to the ITCZ conditions are 
given in the studies of Augstein et al. (1974), Estoque (1975), Godbole and Ghosh (1975), 
and Estoque and Duglas (1978). in Table 5a these data are shown. The relative humidity 
profile in the ITCZ shows a significantly larger value than the climatological average between 
about 800 and 600 mb. In this region of the atmosphere significant deficit in humidity is no- 
ticed when inversion conditions prevail. Utilizing the radiative transfer program we have 
calculated the 9 pm line strength for all these ITCZ cases which are shown in Table 5b. 

The values of 5 and the corresponding w for the ITCZ ase is also plotted in the Fig. 6. 

This Figure clearly shows that the ITCZ C values conform reasonably well with the general 
C vs. w relationship that was derived from mean atmospheric conditions and trade wi..d inver- 
sion models. In all the ITCZ cases the parameter (w - w)/w is negative (see Table 5b) suggest- 
ing that the water vapor content and the depth of the humid layer exceed average conditions. 
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(w — w)/w 


. Figure 11. Relationship between the temperature increase AT, from bottom to top of the inversion and 
the index (w - w)/w based on the data given in Tabic 4. 
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Tabic 5a 

Temperature and relative humidity profiles for observed cases of ITCZ. 
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Table 5b 


Calculated 9 jum line strength, the total water content w, 
and the index (W - w)/W for the ITCZ cases. 



T S (K) 

fi(°C) 

w (g/cm 2 ) 

W (g/cm 2 ) 

(w - w)/W 

JORDAN 

299.15- 

9.92 

4.38 

4.34 

-0.01 

AUGSTE1N 

299.15 


4.94 

4.33 

-0.14 

ESTOQUE a 

297.65 

10.39 

5.26 

4.09 

-0.28 

ESTOQUE b 

298.65 

9.31 

4.58 

4.25 

-0.08 

G0D30LE 

299.65 

11.13 

5.43 

4.41 

-0.23 


We can conclude from the preceding discussion of spectral simulations that the relation- 
ship between £ and w, shown in Fig. 6, is reasonably applicable for the various atmospheric 
conditions we have considered. The parameter (w - w)/w in a general fashion reflects excess 
or deficit of water vapor content in the atmosphere with respect to some mean conditions. 
For positive values this parameter, together with some statistical information, can be related 
to the strength of the trade wind inversion. Negative values of this parameter indicate the 
presence of deep moisture convection such as that produced by ITCZ. 

OBSERVATIONS AND RESULTS 

The Nimbus 4 IRIS obtained spectral measurements for a period of about 9 months from 
April to December 1970. The field of view of IRIS was about 95 km in diameter and this 
instrument had only subsatellitc viewing geometry. As the successive orbits of the satellite 
are about 26° longitude apart and as some data are cloud contaminated, it is not possible to 
construct a global map of the line strength C for each day. It was found necessary to combine 
about 90 days of IRIS data to get a satisfactory global map. For this reason we have divided 
the IRIS data into three periods A pril - June. July - September and October - December, 
and derived three global maps of C. The cloud contamination is eliminated with the help of 
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the brightness temperature in the 1 1 pm window* region. It is assumed that when the window 
temperature exceeds 290 K in the 20°N to 20° S latitude belt over the oceans clouds are ab* 
sent. Similarly a threshold of 285 K and 280 K are applied to the belts 20° - 30° and 30° - 
50°, on either side of the equator, respectively. From the description of the data used in 
this analysis it is clear that some biasing of the information toward clear sky conditions will 
prevail. In addition the composite of data may not lead to a fair spatial and seasonal average. 
Three maps of 5, corresponding to the three time periods mentioned above, derived in this 
fashion, are shown in Figs. 12, 13 and 14. Using *he corresponding 1 1 /am IRIS window data 
and utilizing the split window technique described by Prabhakara et al. (1974) for each one 
of these periods, a sea surface temperature (SST) map is derived. These SST maps are shown 
in Figs. 15, 16 and 17. 

These maps of 8 and SST essentially constitute the basic observational information for the 
subsequent discussion. The data between longitude region of about 100° to 120°E are mis- _ 
sing due to some telemetry limitation . As a consequence the information on this part of the 
Indian Ocean is not contained in our analysis. ) 

t 

The strength of the water vapor lines 8, as shown in the previous section has a quasi-linear 
dependence on the total water. This is not surprising as 8 is a water vapor spectral feature in 
the window region. 

1 

t 

To further validate this result we have used ground truth data obtained from ship and is- 
land radiosonde stations. The satellite 9 pm line strength measurements, within ±1° latitude 
and longitude from the stations, arc used to estimate the total water content according to 
the relationship shown by solid line in Fig. 6. The cloud contamination was avoided by ac- 
cepting only the data in which the surface temperature reported by the ship (or island) dif- 
fered by no more than 1.5°C with respect to satellite derived SST. In this fashion we were 

i 

able to gather 53 cases for comparison, shown in Table 6a and b, covering about 10 month 
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ORIGINAL PAGE IS 
OP POOR QUALITY 


Figure 12. Distribution of the 9pm line strength over the globel oceans (50°N to 40°S) 
derived from the Nimbus 4 IRIS data for the period April, May and June, 1970. 
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Figure 13. Same as Fig. 12 but for the period July, August and September, 1970. 





Figure 14. Same as Fig. 12 but for the period October, November and December, 1970. 





Figure 15. Distribution of the surface temperature over the global oceans f50*N to 40°S) derived from the 1 1pm window region 
observations of the Nimbus 4 IRIS for the period April, May and June, 1970. 
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Figure 16. Same as Fig. 15 but for the period July, August and Septemoer, 1970. 
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Figure 17. Same as Fig. 15 but for the period October, November and December, 1970. 




Table 6a 


Coincident data from ship radiosonde stations and Nimbus ♦ IRIS 


Day 

Ship 


13131 

RK(Po)(%) 

8|RIS(°C) 



Apr. 19 

4 Y M 

275.9 

277.9 

93 

337 

1.16 

0.93 

Apr. 22 

WTK A 

297.0 

299.4 

82 

9.12 

4.11 

3.82 

May? 

4 Y J 

280.4 

281.0 

66 

3.75 

1.32 

0.91 

May 9 

4 YN 

290.5. 

290.9 

68 

4.78 

1.78 

1.90 

May 12 

4 Y M 

278.7 

280.1 

84 

3.03 

1.01 

0.99 

Mav 1 7 

4 YN 

289.9 

291.5 

76 

534 

204 

1.93 

May 21 

■ 

282.4 

283.7 

92 

3.30 

1.15 

1.58 

June 6 

ran 

290.4 

292.9 

63 

5.09 

102 

1.72 

June 12 

mn 

294.3 

292.7 

98 

8.18 

334 

3.46 

June 18 

i 

282.5 

285.1 

76 

3.96 

1.41 

1.35 

June 21 


294.6 

294... 

89 

8.75 

337 

3.71 

June 21 

4 YN 

291.1 

293.1 

82 

5.92 

233 

1.73 

July 2 

4 Y J 

283.9 

284.9 

85 

4.38 

1.61 

1.82 

July 3 

4 Y D 

290.6 

293.1 

92 

539 

2.18 

230 

July 12 

4 Y E 

298.6 

297.7 

77 

8.06 

3.48 

3.50 

July 17 

4 Y V 

295.7 

297.9 

88 

6.44 

2.60 

2.64 

July 31 

4 YD 

293.5 

293.3 

95 

7.12 

2.96 

238 

Aug. 2 

4 YN 

294.8 

296.1 

69 

5.83 

2.28 

133 

Aug. 6 

4 Y B 

280.7 

282.7 

82 

3.17 

1.08 

1.18 

Aug. 7 

4 YC 

283.6 

283.7 

91 

4.87 

1.83 

1.78 

Aug. 9 

4 YN 

295.0 

295.9 

76 

4.96 

1.83 

1.74 

Aug. 1 1 

4 Y D 

294.1 

293.5 

80 

5.96 

2.36 

2.07 

Aug. 23 

4 Y I 

282.4 

285.1 

87 

4.67 

l..« 

1.72 

Ai’.g. 23 

4 YN 

294.7 

294.7 

75 

4.98 

1.88 

1.88 

Aug. 30 

4 YN 

295.6 

295.9 

85 

6.47 

2.62 

2.54 

Aug. 31 

4 YK 

290.8 

293.3 

91 

6.50 

2.63 

3.25 

Sept. 4 

4 YP 

286.0 

283.7 

98 

5.64 

2.20 

2.65 

Sept. 4 

4 Y V 

295.4 

298.7 

80 

8.96 

4.02 

4.14 

Sept. 16 

4 YC 

281.8 

281.9 

81 

4.38 

1.61 

137 

Sept. 27 

4 Y N 

295.1 

294.1 

73 

7.31 

3.06 

2.27 

Oct. 4 

4 Y N 

290.5 

294.9 

. 79 

5.42 

2.08 

2.46 

Oct. 8 

4 Y V 

294.4 

295.3 

65 

5.58 

2.16 

1.50 

Oct. 13 

4 Y J 

286.4 

285.9 

87 

5.26 

2.00 

1.87 

Oct. !8 

4 YN 

294.1 

294.9 

73 

4.90 

1.84 

2.08 

Oct. 2-> 

4 Y P 

279.3 

282.1 

71 

3.18 

108 

0.83 

Oc». 25 

4 Y N 

290.2 

294.3 

73 

5.73 

2.24 

2.12 

Nov. 27 

4 YP 

282.0 

280.9 

80 

3.67 

1.27 

0.91 

Dec. 3 

4 Y V 

291.8 

291.3 

66 

5.00 

1.87 

138 

Dec. 5 

4 Y E 

292.5 

292.2 

72 

6.43 

239 

2.00 

Dec. 26 

4 Y E. 

291.2 

292.7 

85 

6.30 . . j 

233 _ 

__ 2.01 

Jan. 26 

4 Y J 

280.6 

281.5 

75 

3.C7 

1.05 

1 12 











Table 6b 


Same as Table 6a but with radiosonde data from island stations. 


Day 

Station 

m 

m 

RH(p 0 ) 

^IRIS 

W IR1S 

W ISLE 

Apr. 20 

781 18 Turks 

298.3 

296.9 

72 

5.56 

? 16 

2.23 

Apr. 20 

78486 S. Domingo 

_ 

298.4 

296.8 

75 

6.03 

2.41 

2.53 

Apr. 23 

78970 Trinidad 

299.3 

298.1 

90 

7.52 

3.18 

3.57 

Apr. 25 

78526 San Juan 

298.0 

299.3 

67 



3.05 

Apr. 29 

91245 Wake 

298.7 

301.8 

71 

7.38 

3.10 

3.05 

May 3 

78806 Howard 

296.9 

297.7 

93 

10.46 

5.00 

5.07 

May 8 

91165 Lihue 

297.2 

300.2 

72 

6.08 

2.43 

2.57 

May 1 1 

. . J 

781 18 Turks 

297.4 

297.1 

66 

6.14 

2.45 

2.57 

May 13 

91285 Hilo 

299.7 

298.8 

69 

6.16 

2.46 

2.61 

May 18 

78526 San Juan 

298.9 

298.8 

71 

7.30 

3.06 

2.71 

May 18 

78988 Plcsman 

299.5 

300.9 

81 

10.30 

4.90 

5.15 

M-y 21 

78970 Trinidad 

298.7 

300.5 

83 

10.00 

4.70 

4.17 


period over the north Atlantic. The radiosonde data that went into this sample contained at- 
mospheric conditions ranging from inversion to convectively active situations. This compar- 
ison, presented in Fig. 18, reveals that the total water content estimated from the 9 pm line 
strength agrees with the radiosonde measurements within about 1 5%. 

From the relationship, shown in Fig. 6, between £ and w we have constructed three maps. 
Figs. 19, 20, 21, of total water vapor distribution over the global oceans from about 50°N to 
40°S. We are not able to compare this information with a global map of total water derived 
from the conventional data at this time. However Grody et al. (1978) have derived such 
global maps of water vapor using the data from a scanning microwave spectrometer. SCAMS, 














































































W— IRIS (g/cm 2 ) 

Figure 18. Comparison between the water vapor content derived from the 9pm line strength and the 
water vapor content obtained from radiosonde d^ta. 
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Figure 19. Distribution of the total water vapor content (g/cm 2 ) over the global oceans (50°N to 40°S| 
derived from the Nimbus 4 IRIS data over the period April, May and June, 1970, 
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Figure 20. Same 8$ Fig. 19 but for the period July, August and September, 1970. 
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Figure 21. Same at Fig. 19 but for the period October, November and December, 1970. 



that was flown on Nimbus 6. In Fig. 22 one'of these maps' of total water content for the 
period August 18 to September 4, 1975 is shown. We can compare this map with the one 
constructed from IRIS data for the three month period July, August and September, 1970. 
The distribution of the total water vapor in both these maps is very similar. This result sup* 
ports our method for remote sensing of total water vapor over the global oceans. With the 
help of the surface temperature maps and the corresponding maps of the total water content, 
w, we can construct maps of (w - w)/W for the three periods covered by the IRIS observation. 
These maps are shown in Figs. 23, 24, 25. 

In a general fashion these maps reflect the gross structure Oi' the boundary layer over the 
oceans for different seasons. The trade wind inversion associated with the subtropical anti- 
cyclone over the north and south Atlantic and Pacific oceans are clearly delineated by the 
positive value of (w - w)/W. From these positive values we can estimate the increase in tem- 
perature, AT, across the inversion from equation (6). 

Over the Arabian sea strong inversion condition are revealed during the period April - Jure 
as well as during October - December. Apparently the inversion conditions are absent in the 
monsoon period, July - September, giving way to convective activity. 

A vast region of subsidence in the equatorial pacific is revealed with some seasonal changes 
in strength and position. This apparently anomalous phenomenon in the equatorial Pacific 
is known from ship observations (see Riehl, 1954). 

The regions with negative values of (w - w)/w. shown with shading in the figure, indicate 
the presence of eonvectively active areas as explained earlier. Such regions are present where 
the ITC7. generally prevails, as can be seen in the Figs. 23, 24, 25. Convcctivciy active areas 
are also revealed along the course of the Gulf Stream and Kuroshio currents. In the southern 
hemisphere a significantly large belt of eonvectively active zone appears in the October - 
December period. This zone, which is about 10° latitude in width, spans from about 20°S, 
near the east coast of Australia, to about 35°S near Chile in South America. 




Figure 22. Mean precipitable water vapor (g/cm 2 ) derived from Nimbus 6 Scanning Microwave Spectrometer (SCAMS) 
for the period August 18, September 4, 1975. (After Grody, ot al. 1978.) 
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• 23. Distribution of the index, (W - w )/W x 10, over the global oceans derived from the Nimbus 4 IRIS date for 
pril. May and June, 1970. Positive value of this index gives a measure of the temperature increase from bottom tc 
trade wind inversion in °C. Negative values of the index are shown by shading. 
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Figure 24. Seme as Fig. 23 but for the period July, August and September. 1970. 
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Figure 25. Same as Fig. 23 but for the period October, November and December, 1070. 




The significant features of the boundary layer characteristics over - the oceani revealed 
by Figs. 23, 24, 25, are in general agreement with the known climatology. We may conclude 
from this result that it is possible to infer this information from remote measurements of 
total water vapor and surface temperature over the oceans. This inferred information depends 
on one important assumption, that is, the mean water vapor profile over the global oceans 
can be modelled in a simple fashion. Essentially it is assumed that one mean relative humid- 
ity profile, showing a monotonic decrease from surface to 200 mb, represents the average 
condition applicable to all latitudes for different seasons. Departures from such mean con- 
ditions are then associated with the stable condition of inversion or convectively active state. 
For this reason the results obtained in this study are somewhat model dependent. Particularly 
if the mean relative humidity profile, significantly changes as a function of latitude and 
season, the patterns shown in Figs. 23, 24, 25 will be distorted. From the climatological 
data of Nieman it appears that the mean relative humidity profile changes within about 
10n which reflects as a 10% error in w. The error introduced in (w - w)/w from such 
a source has the same magnitude, i.e. 0.1 . Since the range of (W - w )/W shown in the maps 
is much larger than 0.1 we are able to resolve the patterns in a meaningful way. 

The random errors in the remotely sensed SST and w are considerably smoothed in the 
three month averaging process. As a consequence the error introduced in (W - w)/W from 
this source is significantly less than 0.1. 

CONCLUSIONS 

We have demonstrated in this study that it is possible to get a crude information on the 
structure of the boundary layer, over the global oceans, from remote measurements of sur- 
face temperature and total water vapor content. 

Trade wind inversion over the oceans is an integral part of the tropical circulation. The 
studies of Nlalkus (1956), Mak (1976), and Ogura et al. (1977) bring out the importance of_ 



this large scale atmospheric phenomon. The investigation of Reiter (1978) reveals an inti- 
mate connection between the sea surface temperature anomalies and the strength of the trade 
winds. Several studies (see for example Namias, 1978) have revealed the significance of sea 
surface temperature anomalies in producing seasonal climate change. From these studies it 
is clear that the Information we are deriving from remote measurements over the global oceans 
on a seasonal basis, although crude, can be valuable in understanding the ocean atmosphere 
interaction. This knowledge is indispensable in deve’oping seasonal climate prediction 
models. 

The information, obtained with the help of satellite infrared measurements, which we have 
presented here, is some what biased toward clear sky conditions. Microwave remote sensing 
is capable of probing through clouds when they are not precipitating. Thus it appears, when 
the microwave technique to sense the sea surface temperature is developed, the bias intro- 
duced by the clc ads could be reduced. 
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TABLES 


Table 1 Mean profiles of relative humidity and temperature for different latitude belts and 
seasons (from Nieman. 1977). 

Table 2 Mean line strength 5 and total precipitable water w as a function of the surface 
temperature T s for different mean atmospheres. 

Table 3 Three parameters representation of the trade wind inversion (based on Meteor 
data, see Figs. 7 and 8). 

Table 4 Computed 9 pm line strength corresponding to three trade wind inversions (see 
Table 3) for tropics (a), mid-latitude summer (b), and mid-latitude winter (c). 

Table 5 (a) Temperature and relative humidity profiles for observed cases of 1TCZ. 

(b) Calculated 9 pm line strength, the total water content w, and the index (w - 
w)/w for the ITCZ cases. 

Table 6 (a) Coincident data from ship radiosonde stations and Nimbus 4 IRIS. 

(b) Same as Table 6a but with radiosonde data from island stations. 
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FIGURE CAPTIONS 


Figure 1 

Figure 2a. 

Figure 2b. 
Figure 3 

Figure 4 
Figure 5 
Figure 6 

Figure 7 
Figure 8 

Figure 9 


Nimbus 4 IRIS brightness temperature spectrum taken over the equatorial 
Pacific (2.7°N, 140. 1°W, April 20. 1970). Notice the crude structure of the 
water vapor lines at 2.8 cm" 1 resolution in the 9 pm region (1 100- 1230 cm -1 ). 
Water vapor transmission function at 0.1 cm -1 resolution (Kunde and Maguire. 
1974) for 2g/cm 5 of precipitable water. 

Same as Fig. 2a when degraded to 2.8 cm' 1 resolution. 

Mean relative humidity and temperature profiles for tropics and midlatitudes 
summer and winter (from Nieman. 1977). The adopted R.H profile is also 
shown. 

Dependence of the 9 pm mean line strength ? and the total water vapor content 
W on the surface temperature T s over the oceans. 

Weighting functions for the peak at 1 232 cm -1 and valley at 1 225 cm' 1 for the 
three mean atmospheric conditions. 

The relationship between the 9 pin line strength and the total precipitable water 
for mean atmospheric conditions, fot inversion conditions (from Tables 4a, b t 
and c), and 1TCZ cases (sec Tables 5a and b). 

Example of temperature and relative humidity profiles in the trade wind inver- 
sion regime. 

Relationship between the temperature increase. AT, and the decrease in relative 
humidity AR.1I. from bottom to top of inversion, derived from the observation 
of the Meteor expedition (Ticker, 1936). 

Relationship between the height h of the inversion and the increase in temper- 
ature from bottom to top of inversion derived from the observations of the 
Meteor expedition (Fickcr, 1936). 
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Figure 10 

Figure 1 1 

Figure 12 

Figure 13 
Figure 14 
Figure 15 

Figure 16 
Figure 17 
Figure 18 

Figure 19 

Figure 20 
Figure 2 1 
Figure 22 

Figure 23 


Adopted mean relative humidity profile (shown also in Fig. 3) and the models 
of the relative humidity profile for three conditions (see Table 3) of the inversion. 
Relationship between the temperature increase AT, from bottom to top of the 
inversion arid the index iw - w)/tv based on ihe data given in Table 4. 

Distribution of the 9 pm line strength over the global oceans (50°N to 40°S) 
derived from the Nimbus 4 IRIS data for the period April. May and June. 1970. 
Same as Fig. 1 2 but for the period July, August and September, 1970. 

Same as Fig. 12 but for the period October. November and December. 1970. 

Distribution of the surface temperature over the global oceans (50°N to 40°S) 
derived from the 1 1 pm window region observations of the Nimbus 4 IRIS for 
the period April. May and June. 1970. 

Same as Fig. 15 but for the period July. August and September. 1970. 

Same as Fig. 15 but for the period October. November and December, 1970. _ __ 

« 

Comparison between the water vapor content derived from the 9 pm line strength 
and the water vapor content obtained from radiosonde data. 

Distribution of the total water vapor content (g/ctn-) over the global oceans 
(50°N to 40°S) derived from the Nimbus 4 IRIS data over the period April. May 
and June, 1970. 

Same as Fig. 19 but for the period July, August and September, 1970. 

Same as Fig. 19 but for the period October. November and December. 1970. 

Mean precipitablc water vapor (g/cm 1 ) derived from Nimbus 6 Scanning Micro- 
wave Spectrometer (SCAMS) for the period August IS, September 4, 1975. 

(after Grody. et al. l‘)7S.) 

Distribution of the index (vv - w)/w \ 10. over the global oceans derived from 
the Nimbus 4 IRIS data for the period April. May and June. 1970. Positive 
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value of this index gives a measure of the temperature increase from bottom to 
top of trade wind inversio.. .n °C. Negative values of the index are shown by 
shading. 

Figure 24 Same as Fig. 23 but for the period July, August and September, 1970. 

Figure 25 Same as Fig. 23 but for the period October, November and December, 1970. 
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